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Dynamics of spontaneous symmetry breaking and fluctuations in tiie Lipkin-Meshkov-Glick model 
are investigated in a stochastic mean-field approach. Different from the standard mean-field, in the 
stochastic approach, initial state fluctuations, are incorporated. In weak coupling, the approach per- 
fectly reproduces the exact quantal dynamics. On the other hand, for increasing coupling strength, 
above the symmetry breaking threshold, the approach provides description of gross properties (i.e. 
time averaged behavior) of the exact quantal evolution. 
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The mean-field description of a many-body system, i.e. 
the Hartree-Fock (HF) and/or time-dependent Hartree- 
Fock theory (TDHF), provides a simple tool for descrip- 
tions of certain aspects of complex quantum systems [l| . 
For example, in the constrained Hartec-Fock method, by 
explicitly breaking certain symmetries of the underly- 
ing Hamiltonian in static calculation, it is possible to 
describe topology of a quantum phase transitions [^. 
However, it is well-known that the mean-field approxi- 
mation is suitable for the description of mean values of 
one-body observables, while quantum fluctuations of col- 
lective variables are severely underestimated [3| . Numer- 
ous approaches have been proposed either deterministic 
or stochastic to extended mean-field and describe fiuc- 
tuations in collective space 0, Q. Most often, these 
approaches are too complex to be applied in realistic 
situations with actual computational power. A second 
limitation of mean-field dynamics is that it can not de- 
scribe spontaneous symmetry breaking during dynami- 
cal evolution. If certain symmetries are present in the 
initial state, these symmetries are preserved during the 
evolution [l|, 0] . Accordingly, mean- field cannot describe 
physical effects related to spontaneous symmetry break- 
ing including molecule dissociation, spontaneous magne- 
tization, and spontaneous fission in nuclei. 

Both dynamical symmetry breaking and lack of fluctu- 
ations are related to the absence of quantal effects in col- 
lective space and consequently collective motion appears 
nearly classical in the mean-field dynamics. To overcome 
this difficulty, the mean-field approximation should be 
improved by considering a more general wave function by 
coherent superposition of Slater determinants, such as in 
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the time-dependent generator coordinate method [y, 0] . 
However, at present, applications of this method can be 
made only in a very restricted collective space and mostly 
along the adiabatic potential energy surface. Here, we 
employ a stochastic approach, which is much simpler 
than the generator coordinate method, and is based on 
the fact that initial state fluctuations (quantal and ther- 
mal) dominate the fluctuation dynamics at low energies 
[a Q . This idea has been proposed nearly 30 years ago 
by Esbensen et al. in a macroscopic model of nuclear 
reactions [lOl lll|, an d more recently tested in heavy-ion 
fusion reactions [12|- Following a similar idea, recently, 
a stochastic mean-field approach (SMF) has been pro- 
posed [IJ] to treat fluctuations beyond mean-field de- 
scription. In the standard mean-field dynamics, ignor- 
ing quantal and thermal fiuctuations, the initial state is 
specified in a deterministic manner: a given initial state 
leads to a well defined final outcome. In the SMF, on the 
other hand, initial state fluctuations arc incorporated in 
a stochastic approximation. Consequently, an ensemble 
of events are generated starting from a specified distri- 
bution of initial states. It is shown in ref. [l3|, in small 
amplitude limit, that this approach gives rise to the same 
expression for dispersions of one-body observables as the 
one obtained in the variational description of Balian and 
Veneroni (BV) [ij,[i3- In other applications, the average 
version of SMF theory was recently employed jla - llSj to 
successfully reconcile onset of dissipation in TDHF and 
to calculate transport coefficients for relative momentum 
and nucleon-exchange in deep-inelastic heavy-ion colli- 
sions p^ . 

Recently, the variational approach of BV has been 
applied to nuclear reactions [20, l21|. Similarly to the 
standard mean-field description, the approach cannot de- 
scribe spontaneous symmetry breaking mechanism, un- 
less a symmetry-breaking density is used in the varia- 



tional principle. Therefore, the variational approach can 
only provide a poor approximation for dynamical evolu- 
tion in the case of spontaneous symmetry breaking (SSB) 
(see Figure 6 of ref. [22|)- Currently, a realistic descrip- 
tion of spontaneous symmetry breaking in the mean- field 
framework remains an open problem. For this reason it is 
worthwhile to test whether the SMF approach overcomes 
this difficulty. In the SMF approach the initial state is 
not the standard HF state, but specified by a suitable 
distribution. Even if the HF state respects a symmetry, 
in the SMF, this symmetry may be broken initially event 
by event. Consequently one might anticipate that, con- 
trary to the original TDHF and/or BV methods, in the 
SMF approach it may be possible to treat the onset of 
the SSB. We illustrate here that this is indeed the case. 
As a test case, the Lipkin-Meshkov-Glick (LMG) 
Model [3, |23| - [25| is considered here. This model con- 
sists of N particles distributed in two N-fold degenerated 
single-particle states separated by an energy s. The as- 
sociated Hamiltonian is given by (taking h = 1), 



H = e.h - V{Jl - J2), 



(1) 



where V denotes the interaction strength while Ji (i = x, 
y, z), are the quasi-spin operators defined as 
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i(J+ + J_), J, = l(J+-,/_ 
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with J-|_ = X]p=i '^+,p'^-,pi J- = J+ ^'^d where c\_,p and 
c_ are creation operators associated with the upper and 
lower single-particle levels. In the following, energies and 
times are given in e and h/e units respectively. 

This model has the advantage to be exactly solvable 
both in the static [2] and dynamical case [20, [23] . The 
LMG model is known to present a spontaneous symme- 
try breaking in mean-field theory as the interaction V 
increases (see Figure [T] below) . Therefore, this model is 
perfectly suitable to investigate if the SMF approach is 
able to treat the SBB. 

The Hartree-Fock (or Mean-Field) solution is obtained 
by introducing the Slater Determinant trial states written 
as |<f>) — Ilp^ial^p\—) , where the HF single-particle states 
are given by 



-'0,p 



cos(q;)c_ p + sin(a)e*'''c^. p 



vJ 



(3) 



The HF solution is obtained by minimizing the mean- 
field energy with respect to variables a and ip, 



£-iiF[a,ip] 



sN 



|cos(2a) + I sm^{2a) cos(2v3)| (4) 



where x = V{N — l)/e. In Fig. [1] fHF[a, 0] is shown for 
different % parameters. When the strength parameter is 
larger than a critical value (x > 1), the parity symmetry 
is broken in a direction. 
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FIG. 1: (color online). Evolution of the Hartree-Fock energy 
^HF as a function of a for x ~ 0-5 (dashed line), x = 1-8 
(doted line) and x ~ ^ (solid line) for A'^ = 40 particles. 
The arrow indicates the initial condition used in the SMF 
dynamics. 



The mean-field evolution can be formulated either in 



the Schrddinger [27[ or Heisenberg picture. Here, we em- 
ploy the second option. We consider the expectation val- 
ues of the quasi-spin operators ji = {Ji) /N (for i = x, 
y and z). In the mean-field approximation, it is possible 
to derive a set of coupled equations for the expectation 
values of the quasi-spin operators as. 
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Initially, we prepare the system in the state |j, — j), i.e. 
a = 0, which means that all particles are placed in 
the lower single-particle states. This case is indicated 
by an arrow in Figure [T] In this state, initial expecta- 
tion values of quasi-spin components are jz{to) = —1/2, 
jxito) = jy{to) = 0. This initial condition is a station- 
ary solution of Eq. ([5]). When the strength parameter is 
larger that critical value x > 1, the initial state is at the 
saddle point. Since mean-field cannot break the symme- 
try, the system will remain at the saddle point. There- 
fore, it is not possible to describe onset of SSB in the 
standard mean-field framework. This situation is simi- 
lar to the classical object positioned at a = 0. In exact 
quantal description, since the initial state is not an eigen- 
state of the Hamiltonian H , different spin components 
and their correlations change in time. The difference 
between the exact and the mean-field evolution is that 
quantum fluctuations are properly taken into account in 
the exact evolution. 

In the SMF approach, the expectation values of the 
quasi-spin operators obey the same set of equations given 
by Eq. ([5]), except that the initial conditions are differ- 
ent. In order to simulate quantum fluctuations in an ap- 
proximate manner, in the SMF approach [13|, an initial 
ensemble of single-particle density matrices is prepared 
around the same state \j, —j) used in the exact evolu- 
tion. According to the stochastic properties of the initial 



state, it is possible to determine the initial distributions 
of expectation values of quasi-spin operators. We find 
that that the z quasi-spin component is not a fluctuating 
quantity with a mean value jz{to) = ~5- On the other 
hand, the x and y quasi-spin components are uncorre- 
lated Gaussian random numbers with zero mean values, 



j^ito) ^ jiito) 



0, 



and second moments determined by. 



J^ito)j^ito) ^ 3^ito)j^(.to) 



AN' 



(6) 



(7) 



Wc note that even if all trajectories start from the top 
of the energy landscape (the arrow in Fig. [T]) and the 
system has good parity in average, in the SMF evolution 
this symmetry is broken event by event due to non-zero 
values of the spin components along the x and y axis. 

In the SMF, mean values and fluctuations of observ- 
ables are obtained by performing average of expectation 
values over the generated ensemble. The mean values 
of the quasi-spin components and associated dispersions, 
denoted respectively by Ji and Af , are given by 



Mt) = Nj^{t), Afit) = N' (j^r - [jf) . (8) 
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FIG. 2: (color online) Exact evolution of the z quasi-spin 
component obtained when the initial state is \j, —j) for three 
different values of x^ X ~ 0.5 (solid line), x = 1-8 (dotted 
line) and x ~ 5.0 (dashed line) for A'" = 40 particles. The 
corresponding results obtained with the SMF simulations are 
shown with circles, squares and triangles respectively. 

In Fig. [21 the exact and SMF evolutions of the z quasi- 
spin component obtained when the initial state is \j, —j) 
for three different values of %: x — 0-5; X = 1-8 and 
X = 5.0 are shown. In TDHF, for 40 particles, this com- 
ponent remains constant and equal to -20. In both the 
exact results and the SMF simulations, mean values of x 
and y components are zero. In this and following figure, 
the SMF simulation are carried out using a set of 10^ tra- 
jectories. Simulations are performed using Rungc-Kutta 
or order 2 algorithm with a time step of 0.01. 
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FIG. 3: (color online) Exact evolution of dispersions of quasi- 
spin operators obtained when the initial state is \j, —j) for 
three different values of Xi from top to bottom x ~ 0.5 
(top), X = 1-8 (middle) and x ~ 5.0 (bottom) are shown. 
In each case, solid, dashed and dotted lines correspond to 
A^(f), Ay{t) and A?(i), respectively. In each case, results of 
the SMF simulations are shown with triangles (A^), squares 
(Ay) and circles (A^). 



The evolutions of dispersions of quasi-spin components 
obtained in the SMF simulations are shown in Fig. [3l and 
compared with the exact results, Af{t) = (Jf) — {Ji)"^, 
obtained starting from the state |j, — j). We note that, 
in the standard TDHF dynamics, since the state |j, — j), 
does not evolve in time, dispersions of the quasi-spin vari- 
ables remain constant and equal to their initial values, 
A^ = A^ = -J and A^ = 0. Below the critical value of 
the strength parameter (x = 0.5), where energy can be 
regarded as nearly harmonic around a = (see Figure 
[T]), results obtained in the SMF simulations can hardly be 
distinguished from the exact solution. Only a small dif- 
ference is noticeable in the z component. A similar result 
is obtained in ref. [2^ with the BV description. The fact 
that both approaches produce very similar results in the 
harmonic limit is not surprising, since it was shown that 
they both contain the same physical information in this 
limit [1J|. Above the critical strength (x = 1.8 and 5) 
(middle and bottom panel of Figure [3|) , the B V descrip- 



tion has been shown to lead to very bad results |22| , when 
calculations start from the same initial condition. Here, 
we see that the SMF approach provides a fairly good re- 
production of gross properties of the exact dynamics. In 
particular, during the early times, the SMF simulations 
can not be distinguished from the exact evolution. Dur- 
ing the long time evolution, the SMF simulations describe 
time-averaged behavior of the exact dynamics very well. 
As seen from Fig. [2j a similar agreement is obtained for 
the mean value of the z component of quasi-spin for all 
values of the strength parameter x- The present example, 
clearly demonstrates the ability of the SMF approach to 
describe gross properties of mean values and fluctuations 
for any strength of the interaction. 

It is well known that the standard mean-field theory 
provides a poor description for fiuctuations of collective 
motion, and it essentially treats the collective motion in 
a classical approximation. The SMF approach makes 
an attempt to correct this shortcoming by incorporat- 
ing quantal and thermal fluctuations in the initial state. 
In this work, we test the approach in the LMG model. 
As seen in Fig. [5] and Fig. [31 the SMF simulations 
provide nearly perfect description for non-trivial oscil- 
lations during early evolution of mean values and dis- 
persions of quasi-spin operators. Over the long time in- 



terval, simulations also provide a satisfactory description 
for the gross properties, i.e., time averaged behavior of 
the mean values and dispersions of quasi-spin operators. 
We should note that we do not expect that such a simple 
SMF approach provides a detailed quantum mechanical 
feature of the evolution. In particular, possible interfer- 
ences between different trajectories are neglected as well 
as possible tunneling effects. Nevertheless, the stochastic 
method presented here provides a suitable framework be- 
yond mean-field for describing dynamics of fluctuations 
and for understanding spontaneous symmetry breaking 
in complex quantum systems from a quasi-classical per- 
spective. 
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